Imperfect measurements settings: 
implications on quantum state tomography and entanglement witnesses 
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Reliable and well-characterized quantum resources are indispensable ingredients in quantum in- 
formation processing. Typically, in a realistic characterization of these resources, apparatuses come 
with intrinsic uncertainties that can manifest themselves in the form of systematic errors. While 
systematic errors are generally accounted for through careful calibration, the effect of remaining 
imperfections on the characterization of quantum resources has been largely overlooked in the liter- 
ature. In this paper, we investigate the effect of systematic errors that arise from imperfect alignment 
of measurement bases — an error that can conceivably take place due to the limited controlability 
of measurement devices. We show that characterization of quantum resources using quantum state 
tomography or entanglement witnesses can be undermined with an amount of such imprecision that 
is not uncommon in laboratories. Curiously, for quantum state tomography, we find that having 
entanglement can help to reduce the susceptibility to this kind of error. We also briefly discuss how 
a given entanglement witness can be modified to incorporate the effect of such errors. 



I. INTRODUCTION 

The advent of quantum information science has 
brought inspiring opportunities for information process- 
ing [1(. At the heart of all these information process- 
ing protocols is the encoding of specific information in 
quantum systems and the ability to perform some specific 
measurements — these features arc notably important, 
e.g., in measurement based quantum computation [2Jj. 
Evidently, real life preparation of specific quantum states 
is never ideal and thus their reliable characterization is 
crucial for the implementation of these protocols. 

For a complete characterization of quantum state, one 
uses the technique of quantum state tomography (see, 
e.g. Refs. 0, |j]) whereas for the purpose of entangle- 
ment verification, the technique of measuring entangle- 
ment witnesses [5|48| is widely employed. A common fea- 
ture of these techniques is that they require measure- 
ments to be carried in a number of different settings. 
While these settings can be theoretically established eas- 
ily, their experimental implementation may differ from 
the theoretical prescription, or come with intrinsic un- 
certainty, thus contributing to non-negligible systematic 
error. For instance, the measurement on a polarization 
qubit cannot be more precise than the intrinsic uncer- 
tainty of the polarization rotator used (typically of the 
order of 1° ~ 2° in real spacqB). Likewise, the precision 
of measurements on two-level atoms is limited by the ef- 
fective phase and intensity uncertainties of the laser pulse 
experienced by the atom. 

In this regard, it came as a surprise that the intrinsic 
uncertainty or systematic error present in measurement 
devices is — to our knowledge — hardly reported in ex- 
perimental findings, clearly in stark contrast with statis- 
tical error j9|. Also, the implication of imperfect devices 



1 This translates to an uncertainty of 2° 
Bloch sphere. 
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seems hardly investigated bey ond a relatively small num- 
ber of research articles [3, Hfl, [ll|] . Of course, with careful 
calibration, systematic errors can usually be detected and 
reduced, see, e.g., Refs. |12h16J . However, it is important 
to note that even after careful calibration, measurement 
devices are after all never perfect and the intrinsic un- 
certainties can still manifest themselves in the form of 
bounded systematic errors. 

The main purpose of this paper is to present con- 
crete evidence showing that the potential implications 
of overlooking systematic errors can be significant in 
the characterization of quantum states. We illustrate 
this by considering a specific kind of systematic error 
that arises from misaligned measurement bases, and il- 
lustrating its effect on two commonly employed meth- 
ods for characterizing quantum states, namely, quantum 
state tomography and the evaluation of entanglement 
witnesses [8J]. Our analysis therefore complements the 
approach of Ref. [17| > which allows the detection of sys- 
tematic error from experimental data. 

The paper is structured as follows. We begin in Sec. |TT] 
by explaining the misalignment systematic error that 
arises from imperfect measurement settings. The nota- 
tions that we are going to use in the text will be in- 
troduced in the same section. In Sec. IIII1 we present 
the effect of such systematic errors on quantum state to- 
mography, in particular the fidelity of the reconstructed 
state with respect to the actual state prepared. Next, in 
Sec. IIV1 we illustrate the effect of misalignment error on 
the evaluation of a family of genuine multipartite entan- 
glement witnesses. In both Sec. IHII and Sec. IIV1 we also 
discuss how these effects can be compensated when the 
amount of systematic error (uncertainty) is known. We 
conclude in Sec. [V] with a summary of main results and 
some possibilities for future research. Technical details 
related to the main results can be found in the Appen- 
dices. 



II. MISALIGNED BASES FROM 
IMPERFECT MEASUREMENTS 

To study the effect of misalignment systematic errors 
on the characterization of quantum states, we consider 
the typical scenario where n spatially separated qubits 
can each be measured locally in a number of different 
bases (settings). For the benefit of subsequent discussion, 
we remind that any of these local measurements can be 
described in terms of a 3-dimensional unit vector on the 
Bloch sphere. More explicitly, we shall denote by 
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the fc-th qubit observable! to be measured on the j-th 
qubit and mi the corresponding Bloch vector. A mis- 
alignment error can then be defined as follows. 



Definition 1. A misalignment error is a systematic er- 
ror that arises from imperfect measurement settings, i.e., 
the actual observable measured N^ = n)! -cr differs from 
the intended measurement setting Mf, for at least some 
value of j and k. 

Henceforth, we shall quantify the amount of error 
present by the quantity 



- (?) ~(i) 
maxacos I m^ ■ nj: 

3>k 
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Geometrically, this means that among all the measure- 
ment settings chosen by all the n parties, the maximal 
angular deviation^ of the actual measurement directions 
from the intended ones are at most e. 

To simplify the subsequent discussion, we shall also 
assume that all outcome probabilities can be estimated 
with negligible statistical error and that the actual mea- 
surement bases can always be described by fif! in all 
runs of the experiments. Next, we look into the effect 
of such misalignment error on some commonly employed 
protocols used in the characterization of a quantum state 
(see also Figure Q]). 
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FIG. 1. A schematic representation of the procedures in- 
volved in the characterization of a quantum state and how 
misalignment error can affect the resulting characterization. 
(1) The characterization protocol specifies measurement in a 
certain bases {intended measurements). (2) Due to imperfect 
measurement devices, misalignment errors creeps in during 
the experiments and the actual measurements performed dif- 
fer from the intended ones. (3) The resulting measurement 
statistics are used to compute a given figure of merit, such as 
the fidelity or the expectation value of an entanglement wit- 
ness. (4) Imprecision of the measurement e then translates 
into additional uncertainty in the final figure of merit. 



Throughout, we shall assume that the qubit tomogra- 
phy is intended to be carried out in the standard Pauli 
bases, i.e., 
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for all parties. Moreover, we shall quantify the effect of 
misalignment errors on quantum state tomography using 
the Uhlmann-Jozsa fidelity [Ija, [lj| between the actual 
state t and the reconstructed state p, i.e., 



F(t,p) = f tryVrp\/rJ 



(4) 



It is worth noting that when cither p or r is a pure state, 
the expression above reduces to J-{t, p) = tr(pr). 

Clearly, the smaller is the value of F(t,p), the more 
drastic is the effect of misalignment error on quantum 
state tomography. 



A. Single qubit state tomography 



III. IMPLICATIONS ON QUANTUM STATE 
TOMOGRAPHY 



Quantum state tomography is the process in which 
many copies of a quantum state are measured in a set of 
tomographically complete bases, followed by some state 
reconstruction algorithm using the measurement statis- 
tics and the presumed knowledge of the measurement 
bases [3|, Q, LLI| • In this section, we illustrate the effect 
of misalignment systematic error on the tomography of 
n-partite qubit states. 



2 Here, a = {a x ,(Ty,a z ) denotes the vector of Pauli matrices. 

3 As measured in the Bloch sphere. 



Let us begin with the simplest example of a single- 
qubit state tomography. The pedagogical example given 
below will also serve to remind the key features involved 
in some of the standard state reconstruction techniques, 
such as linear inversion and maximum-likelihood estima- 
tion. 



1. A simple example of erroneous state reconstruction 
starting from a pure state 

Consider a source that produces a quantum state r as 
parametrized by the Bloch vector t: 



1 + t- a 



(5) 



Suppose now that a qubit state tomography is to be car- 
ried out for this source via the intended measurement^ 
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Mi = a x , M 2 = (T y , M 3 = cr z , 



(6) 



whereas in reality, due to the presence of misalignment 
errors, the actual observables measured arc described in- 
stead by {N k }. 

From Born's rule, we can compute the outcome prob- 
ability for the k-th measurement setting as: 

P(±l|fc)=tr(r 1± iVfc ), for k =1,2, 3. (7) 



The essence of state reconstruction is to find a legitimate 
density matrix p, referred as the reconstructed state such 
that 



P(±l|fc) =tr [ p 



1± M fc 



(8) 



Since this amounts to solving a set of equations that are 
linear in the measurement statistics, this procedure of 
solving for the reconstructed state p is also known as lin- 
ear inversion. Note that the reconstruction is done using 
the ideal description given in Eq. ^. Evidently, if the 
actual misalignment error was detected, the reconstruc- 
tion procedure could be corrected by replacing {M k } with 
{N k j in Eq. ©. 

Specifically, imagine that the actual state prepared is 
r = |V , ^)('0^l: where |^) is the positive eigenstate of 
s ■ a with s 



^(W) T = t: 
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and the actual tomography measurement directions di- 
verge uniformly from the intended directions as e in- 
creases (see Figure [2]), i.e., 



n\ = 



where 
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c e = cose, s e = sine. 
It now follows from Eq. ©-Eq. flTU]) that 



(10) 



(11) 




FIG. 2. Intended and actual measurement directions for 
the tomography of \ip s )- The actual measurement direc- 
tions {rife} given in Eq. (|10[) correspond to the triad of {rhk} 
"opened" uniformly towards — s — — 7= (1, 1, 1) while satisfy- 
ing acos(m fc -h k ) — e. 



The Bloch vector t = s of r = \"4>f)(ipf\ is thus an eigen- 
vector of the above linear transformation, with eigenvalue 
c e — \/~2s e . Hence, as long as \c E — y/2s e \ < 1, which takes 
place for < e < acos| w 70°, the reconstructed state 
obtained by solving Eq. (|8]) is always a legitimate quan- 
tum state. It is then straightforward to verify that p can 
be written as a convex mixture: 



p = f\riM\ + (i-f)\^)&i 



where 



/(e) = ^(r,p) = -(l 
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(13) 



(14) 



is simply the fidelity of the reconstructed state p with 
respect to the actual state r = \ipf}(fp^\ [cf. Eq. (Q}]. 

This implies, for instance, that with a 2° misalignment 
in rhk but everything else perfect, the erroneously recon- 
structed state still only has 97.5% fidelity with respect 
to the actual state (Figure [3]). In this regard, note that 
Eq. (fT4|) actually represents the worst-case fidelity for 
any actual state r that is pure and where the intended 
tomographic measurements are given by Eq. (j6|). The 
proof of this is somewhat involved and is relegated to Ap- 
pendix [XTJ More generally, to study the effect on fidelity 
for small e, we shall introduce the notion of susceptibility 
to misalignment errors, defined as: 
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In particular <S 7^ shows that misalignment errors have 
a first-order effect on the fidelity. 
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(c e - V2s e )s. (12) 



For simplicity, we omit all superscripts in the single-qubit sce- 



2. Effect on state reconstruction starting from a general 
qubit state 

Obviously, in a realistic experimental situation, we do 
not expect any source to produce a pure qubit state. A 
relevant problem to determine is thus whether the ex- 
pression given in Eq. (fT4")) still represents the worst-case 
fidelity even in this more general scenario. 
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FIG. 3. Fidelity of the reconstructed state with respect to 
the initial state \ipt)® n when the actual measurement settings 
are those described in Figure [2] and the intended measurement 
bases are those given in Eq. @. The vertical dotted line 
corresponds to e = 2° as discussed in the text. 



Before answering the above question, it is important to 
understand that the linear inversion process described in 
Sec. lIII ATl does not always lead to a physical state p. For 
instance, in the example given above, if e > acos^, the 
reconstructed Bloch vector as given by Eq. (fT2j) would 
have length greater than 1, and thus corresponds to an 
unphysical state p. 

A commonly employed technique to circumvent this 
kind of problem is to make use of the maximum- 
likelihood estimation (MLE) technique introduced by 
Hradil [201 ■ From the measurement statistics and the 
supposed knowledge of the measurement bases, this tech- 
nique seeks to find a physical quantum state that maxi- 
mizes the (log) likelihood function, and hence determines 
the quantum state that is most likely to give rise to the 
experimental data. 

For a general mixed qubit state and a set of 3 mea- 
surement directions {fik} satisfying Eq. ([2]), one does 
not always obtain a physical state via linear inversion. 
Nonetheless, we prove in Appendix I A II that with MLE, 
the fidelity of the reconstructed state cannot be worse 
than that given in Eq. (fl"4)) . i.e., 



T(r',p')>f(s), 



(16) 



for any qubit state t' and the corresponding state p' re- 
constructed from the MLE algorithm. 



B. Multiqubit state tomography 

Let us now study the effect of misalignment error on 
the tomography of multiqubit states. As we will see be- 
low, in the reconstruction of multiqubit states, entangle- 
ment also plays a nontrivial role in combating the effect 
of misalignment errors. 



1. Multiqubit product states 

To start off, note that the simple one-qubit example 
given above can be easily generalized to the n-partite 
scenario if the source actually produces a product state, 
i.e., 
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(17) 



j'=i 



To see this, we note that the product nature of quan- 
tum states is preserved by the MLE state reconstruction 
technique (for a proof of this, see, e.g., Appendix I A 2[) . 
Moreover, the reconstruction procedure can be carried 
out independently for each qubit. Thus, if we define anal- 
ogously f n {e) the worse-case fidelity in the n-partite case, 
it follows that 
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(19) 



Note that the inequality in Eq. (JT5)) is saturated by con- 
sidering r = I^X^I®™ and where all its constituents 
are measured with axes defined in Eq. (|10l) (see also Fig- 
urc[2|). Clearly, this shows that the effect of misalignment 
error may accumulate with the number of parties. The 
fidelity itself /„ (e) as a function of e for n < 4 is plotted 
in Figure [3] 



2. Two- qubit- entangled states 

Evidently, in the context of quantum information pro- 
cessing, it is arguably more relevant to look into the ro- 
bustness of entangled states with respect to the aforemen- 
tioned systematic errors. To this end, we have performed 
numerical optimization to determine — for small s and 
for fixed amount of entanglement (as parametrized by 
a E [0, j J) — the worst-case fidelity F(\ip a ),p) by vary- 
ing over the misaligned measurement settings (as param- 
eterized by n£ ) and arbitrary qubit basis states |^- ) 
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LQ#1 )l^2 



(20) 



In our optimization^ we assume Eq. ([3]) and focus on 
small error, namely, e < ^Sn ^° determine the worst- 
case fidelity .F(|-0 Q ), p) numerically for e in this domain. 



5 To make the optimization more efficient and robust, we have also 
provided the gradient of the objection function (with respect to 
the parameters of the problem) to the optimization solver. 
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We then estimate numerically the susceptibility, i.e., the 
initial slope S (a) = \7 £ J : (\ip a ) 7 p)\ E= o- The results are 
shown in Figure. |4] Interestingly, our results show that 
in the worst-case scenario, pure product states are the 
least robust against systematic errors that arise from mis- 
aligned measurements. In the two-qubit case (n = 2), 
Eq. (|18l) thus provides the worst-case rate of decrease of 
fidelity with respect to e for small e. It is also interesting 
to note that maximally entangled two-qubit pure state 
does not appear to be the most robust against this kind 
of error. 
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FIG. 4. Susceptibility of 2-qubit pure state \ip a ) to misalign- 
ment error e as a function of the entan glem ent present in 
\ip a ) (parameterized by the concurrence [211 ] of \ip a ))- For 
given a, S(a) gives the rate of decrease of the fidelity with 
respect to e as e — ► 0; S(0) is the corresponding intial slope 
for pure product state. Fhe 257 numerical data points ob- 
tained from 4 x 10 optimizations are plotted in a solid line. 
Fhe dashed lines represent segments of the plot that can be 
very well approximated using the explicit parameterizations 
given in Appendix I A 31 

What gives entangled state more resistance to this kind 
of systematic error in the worst-case scenario? Our intu- 
ition is that uncorrelated, local misalignment errors have 
mostly local effect. Here, the misalignment errors are 
uncorrelated in the sense that for each intended mcasurc- 

(2) 

ment direction, say, for the second party rh k , its actual, 

deviated measurement direction n. is independent of 
the choice of measurement k\ by the first party. To verify 
this intuition, we have performed similar analysis allow- 
ing the actual measurement direction to vary depending 
on the choice of measurement of the other party. Indeed, 
it turns out that pure product state is no longer the most 
fragile one against misalignment error in this more gen- 
eral scenario. More details on this analysis can be found 
in Appendix I A 41 

Coming back to the uncorrelated case, we note that 
for small amount of entanglement, say, a < |^, the re- 
constructed state p loses its fidelity with respect to the 
actual state \ipa) — as quantified by 1— J-(\ip a ), p) — pre- 
dominantly via terms that are proportional to the length 
of the Bloch vector of the reduced density matrix. Since 
this length cos la shrinks as a increases, clearly, among 
all the weakly entangled two-qubit pure states, the pure 
product state is the most susceptible to this kind of sys- 



tematic error. A more formal analysis of this is given in 
Appendix IA3I . 



IV. 



IMPLICATIONS ON ENTANGLEMENT 
CERTIFICATION 



While standard quantum state tomography can be car- 
ried out for a system involving a small number of qubits, 
in the realms where quantum information processing is 
advantageous against its classical counterpart, this com- 
plete characterization is practically infeasible (see, how- 
ever, Ref. |22()- Next, we shall look at the implication of 
misaligned measurements on partial characterization of 
quantum state via entanglement witness. 



A. Preliminaries 

A witness W for genuine ri-partite entanglement is a 
Hermitian observable that satisfies 



tr (Wp h 



iscp. 



> 0, (21) 

for all biseparable states p(bi-)sep. but is violated by at 
least some (genuinely n-partite) entangled states [3, |8[ . 
An optimized entanglement witness W, moreover, satis- 
fies the property that there must exist biseparable quan- 
tum state Pbisep. such that the defining inequality, cf. 
Eq. (|2"Tj) is saturated. Geometrically, this means that the 
separating hyperplane defined by W is actually tangen- 
tial to the set of biseparable states (see Figure [5]). 
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FIG. 5. A schematic diagram showing the effect of misalign- 
ment error on the evaluation of an entanglement witness. An 
optimized entanglement witness W is tangent to the bound- 
ary of the set of biseparable states, (a) When evaluated using 
misaligned measurements, the witness can cross the bound- 
ary (as W e ) and thus some biseparable states may appear to 
be genuinely n-partite entangled, (b) To correct the prob- 
lem, one can evaluate the correction factor w(e) such that no 
biseparable state give a false positive result. 

Though being more economical in terms of resource 
requirements, we shall demonstrate below that entangle- 
ment certification via (optimized) entanglement witness 
is relatively fragile against misalignment systematic er- 
rors (see also Figure [S]). For definiteness, we assume 
in subsequent analysis that a linear entanglement wit- 
ness W is intended to be evaluated by measuring local 



observables M. 
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a. And as with the rest of - along the following directions on the Bloch sphere 



the paper, we assume that due to the presence of mis- 
alignment errors, the actual local observable measured is 
a, where the angular deviation of fif! from 

m]^' is bounded by e, cf. Eq. ©. 

To incorporate the effect of uncorrelated misalignment 
error, one can first determine the correction factor 



Ai) 



w(e) = min min tr (W e 

W e Pbisep. 



Pbiscp. ) 



(22) 



where the minimization of W e is to be carried out over 
all possible Hermitian observables W e satisfying the con- 
straint given in Eq. ([2]). We write W 6 * the Hermitian 
observable giving the minimal value of u;(e)|j 

The function w(e) thus gives the worst-case value of the 
witness W with respect to all biscparable states in the 
presence of bounded misalignment error e. If e is known, 
the witness W can then be modified in the following way 



such that 



W^W' = W-wie)^ 



tr (WVbisep.) > 0, 



(23) 

(24) 
even if we allow 



holds true for all biseparable states pbisop 
misalignment error bounded by e, sec Figure ([5]). 



B. A bipartite entanglement witness and its 

correction factor in the presence of bounded 

misalignment error 

Let us now look at some explicit examples. Con- 
sider the following two-qubit entanglement witness con- 
structed from the singlet state I* - ), 

w*- = \i® 2 -|*-)(*- 1 = \^ 2 + \ E °*®^ 

A;— x,y,z 

(25) 
where 1 is the 2x2 identity matrix. Clearly, a natu- 
ral way to evaluate this witness experimentally involves 

(?) 
measurements in the Pauli bases, i.e., with M? given 

by Eq. ©. 

Imagine now a physical system prepared in the sepa- 
rable state: 



|V;> = cos 2 x (|0) + e** tan X |l» (tan X |0) + e~^|l)) , 

(26) 

where x = as °2 an0 -) histead of the Pauli bases, mea- 
surements were made — due to misaligned measurements 



Note that by convexity of the set of biseparable states, it suf- 
fices to consider pure biseparable quantum state l^bisep) in the 
minimization of Eq. (122 I t . 
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An intended measurement on W^- using \ip) therefore 
results in the measurement of 



fe=l 



<<$ 



^nP-a®nW-a, 



which gives an expectation value of 

(^|W|_ \ip) = I (cos 2e - 2\/2sin 2e - 1 
8 V 



(28) 



(29) 



which is negative for all < e < ?. In other words, as 
soon as e > 0, an evaluation of the above entanglement 
witness W^,- using the separable state \ip) will always 
give an affirmative, but erroneous certification that the 
state is entangled. Numerically, the above strategy also 
corresponds to the minimal value that we have found 
for the optimization specified in Eq. (|2"2"T) . Therefore, for 
the witness W*- , our result suggests that the correction 
factor w(e) is given by Eq. (|29l) . 



C. A witness for genuine n-qubit entanglement and 

its correction factor in the presence of bounded 

misalignment error 

Likewise, for the detection of genuine multipartite en- 
tanglement, let us consider the following n-partite entan- 
glement witness |23j : 
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Wghz = -1®" - |GHZXGHZ|, 



(30) 



where I GHZ) 



775(10) 
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l)® n is the n-partite 



Greenberger-Horne-Zeilinger state. An economical way 
to measure this witness is to have all the n parties per- 
forming the same measurements j24j . i.e., 
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and M, 
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cos — G x 
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kn 
sin — a. 
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— <7 y , for k = 1, . . .,n, (31) 
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a,. The measurement statistics on these 



settings can then be combined to give the desired expec- 
tation value of Wghz hi the following way: 
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(32) 
In what follows, we provide estimates of the correction 
factor wghz(£) of W|ghz) — as determined by numerical 
optimization — for n > 4, separating the cases of n even 
and n odd. This correction factor u>ghz(e) is plotted for 
n = 3, . . . , 8 in Figure [f)J 
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FIG. 6. The worst expectation value of W E found using 
bispearable states and systematic error e < 10°. When e = 0, 
W e reduces to W|ghz) given in Eq. (|32p . The curves for n > 4 
are computed using Eqs. (J33J and (|37[) . whereas the curve for 
n = 3 has been obtained numerically. 




FIG. 7. Evaluating the entanglement witness Wghz with the 
intended measurement settings M^ on the XY plane (drawn 
here for n = 4 parties). These settings correspond to equally 
spaced measurement directions m); on the XY plane, with 
neighboring ones separated by an angle — . The systematic 
errors considered here correspond to having these rays moving 
in pairs towards each other as e increases. Measurement di- 
rections are deviated according to the black arrows for qubits 
numbered j — 1,2 and to the grey arrows for j = 3, 4. Mea- 
surement in the <r z basis is assumed to stay unperturbed. 



1. Estimated correction factor for even n > 4 

For W|ghz) with even n > 4 and e < £-, numerical 
optimizations suggest that the correction factor is given 
by: 



n oven/ 
W GHZ I 



■ sin ne. 



(33) 



This can be achieved by considering the n-qubit bisepa- 
rable state 

|V>„) = I(|o>®*+e**|l)®* )®(|0}®§+e- i 5|l)®f), (34) 

for even n > 4, with the following misaligned observables: 



ti<- 



cos 



N 



and 



Nl j> * } = cos 



kn , „ i, 
— + (-l) k e 




a x + sin 


n 





kit , . fc 
— + (-l) k e 
n 



(35) 



^ - (-l) fe e 
n 



a r + sin 



--(-l) fc e 



(36) 

(see Figure [7]) . Clearly, Eq. (|33|) is negative as soon as 
e > 0. Thus, as with the two-qubit example given above, 
measuring the biseparable state \ip n ) for these imper- 
fectly implemented witnesses always results in an erro- 
neous certification of the non-separability of the state. 



2. Estimated correction factor for odd n > 5 

For W|ghz) w ith odd n > 5 and e < ^, numerical 
optimizations suggest that the correction factor is given 



by: 



«>) = Tn 



n — 2 — (n — 1) cos e + cos ne • 



sin ne 
tan ■£- 

2n 



(37) 
This can be achieved using the biseparable pure state 

|V>„) = |iM®rV»+)> 



l^_> = -L(|o>®»-+. 
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IV>+> = -3 ( |0>' 



§n+ _|_ p^' (3 " 4 t,' /) ' r |1\®" 



|i\ig)n_ 

i \0Sm+ 



(38) 
(39) 



where n± = ^^- and u = (— 1)™- . We choose the follow- 
ing local measurement settings for k = 1, . . . , n: 



N^ = cos 



kir (j\ 
— +9i°'e 
n 



Ot + sin 



kn 



or 



— + 9k £ 



<Ty, (40) 



keeping N^ = M^_ x = a z and defining #£" 
^(-l) fc sign[(fc-n + )(j-f)]. 



V. CONCLUSION 

Intrinsic uncertainties in measurement devices, which 
can manifest themselves in the form of misalignment sys- 
tematic errors represent an unavoidable part of any real- 
life quantum experiment. In this paper, we show by 
explicit examples that the procedure of characterizing 
quantum resources using state tomography or entangle- 
ment witnesses can be considerably affected when such 
systematic errors are not properly taken care of. For 
example, when considering pure two-qubit state tomog- 
raphy, every single degree of misalignment on the Bloch 
sphere can potentially lead to ~ 1% decrease in the fi- 
delity of the reconstructed state. For general product 



state, the worst loss of fidelity has a scaling that is linear 
in n, making them increasingly sensitive to such system- 
atic error. 

For the verification of entanglement via an optimized 
entanglement witness, we show that an erroneous certifi- 
cation could arise whenever there is nonzero misalign- 
ment error. While our demonstration was made for 
specific entanglement witnesses, it should be empha- 
sized that all non-device-independent entanglement wit- 
nesses [25| are potentially susceptible to this kind of im- 
perfection, and the procedure we followed can also be 
applied to them. But all is not lost, the effect of mis- 
alignment error, as we have demonstrated, can be in- 
corporated by modifying a given entanglement witness. 
In this regard, it could be interesting to understand 
the amount of potential misalignment systematic error 
present in some of the state-of-the-art characterizations 
of quantum resources, such as those in Refs. [2a - [28l |. 
Alternatively, entanglement verification can also be car- 
ried out without such characterization by implemen ting 
device-independent entanglement witnesses [25l . l29l . l30j 
provided by Bell-like inequalities. 

Let us now comment on some further possibilities for 
future research. Firstly, for quantum state tomography, 
our analysis of entangled two-qubit states focused on 
pure states; based on the numerical optimizations that 
we have done, we conjecture that the bound given in 
Eq. (18) holds for mixed two-qubit states as well. Ob- 
viously, similar studies for n-qubit systems should be 
carried out for n 3. The increased resistance that we 
have observed in entangled two-qubit states against mis- 
alignment error suggests that entanglement also plays a 
nontrivial role in quantum state estimation something 
that deserves to be understood better. For instance, 
it would be interesting to see what role entanglement 
plays when considering other imperfections, such as mis- 
matched detector efficiencies. We remind also that our 
analysis on the effect of misaligned measurements was 
carried out at an abstract level where, in particular, each 
measurement basis can be misaligned differently but in 
an uncorrelated manner (see, however, Appendix A4). 
In practice, typical errors present in particular experi- 
mental setups may be more/less general than considered 
here, leading to larger/smaller effects. It would thus be 
interesting to adapt the analysis that we have presented 
here for some actual physical system (e.g. superconduct- 
ing qubits, qubits in ion traps) and see how the effect 
changes. 

Clearly, our work only marks the beginning of a deeper 
understanding how imperfect devices can affect real-life 
characterization of quantum resources. In the long run, it 
is clearly desirable to develop a general method for com- 
puting the additional uncertainty that should be incor- 
porated in any figure of merit as a result of any given im- 
precision in the measurement device. The joint effect of 
imperfect devices and finite statistical error is evidently 
also a relevant question that needs to be addressed. 

Of course, it is also of general interest to understand 



how imperfect measurement settings directly affect quan- 
tum information processing tasks, which evidently re- 
quire more than well-characterized quantum resource. To 
this end, we note that the effect of imperfect measure- 
ment settings on measure and prepare quantum key dis- 
tribution protocol is investigated in parallel in Ref. [3l| . 
Acknowledgments. Wc acknowledge useful discus- 
sions with Stcfano Pironio and Clara Osorio. This work 
was supported by the Swiss NCCRs QP and QSIT, the 
CHIST-ERA DIQIP and the European ERC-AG QORE. 



Appendix A: Miscellaneous details related to 
quantum state tomography 

1. Tomography of a single qubit mixed state 

Under the assumption that the misalignment system- 
atic error is upper-bounded by e, we show in this Ap- 
pendix that the minimal fidelity of the reconstructed 
single-qubit state p with respect to the actual state r, 
i.e., F(r,p) is indeed given by Eq. (fblj) . Throughout 

this Appendix, we assume that e < acosw| « 35° and 
that the state is reconstructed via linear inversion when- 
ever possible, or otherwise via the maximum-likelihood 
(MLE) estimation technique. 



a. Outcome data compatible with r, e 

We start by noting that for the purpose of state recon- 
struction, instead of the outcome probability distribution 
computed in Eq. (J7]), wc can just as well work with the 
vector c £ R 3 defined as follows Q 

c such that c k = P{+l\k) - P(-l\k) = tr (rN k ). 

(Al) 
In the absence of misalignment error, i.e., when N k = M k 
for all fc, c is simply the vector of average values with 
respect to the Pauli matrices. 

We now characterize the set C of c obtainable from the 
actual state r = (1 + 1 ■ <?)/2 and bounded misalignment 
error e. Using Eqs. (JSJ and (|A1[) . wc see that c and t are 
related through a linear transformation: 



tr [(hi 




(A2) 



where $ is a real matrix that encodes the actual mea- 
surement directions, h, x is the x-component of the unit 



7 There is a one-to-one correspondence between the components 
of c and the measured probability distribution P(±l|fc). 



vector hk\ njF and n fc arc analogously defined. To 
parametrize C(r,e), we recall from Eq. @ that the mis- 
alignment errors are bounded such that n^ satisfies 



cose < rife • rhk- 



(A3) 



Additionally, we observe from Eq. (|A2[) that each com- 
ponent of the vector c can be written as Cfc = hk-t. Thus 
for given t, each Ck is constrained with an interval. More 
precisely, the set C(r,e) is a box whose boundaries are 
specified by vectors saturating inequality (|A3j) . 

To simplify the computation, we now show that C(t, e) 
is contained inside a ball B(t, tX) of radius tX centered at 
T(see Figure E]), i.e., B(t,tX) D C(r,e), where 



A = 1 — cose + v2sine, 
and e < 35° ensures < A < 1. 



Block sphere B(t, tX) 



(A4) 





FIG. 9. Reconstructed state r * using MLE, compatible with 
outcome data c outside the Bloch sphere. We show in the 
text that cis contained inside a ball of radius tX centered at t 
(dashed line). The reconstructed Bloch vector r* lies on the 
surface of the Bloch sphere. Moreover, since the function Of) 
is strictly concave — schematic representation of the contour 
lines of £(r) are plotted using dots — r * is constrained on a 
spherical cap such that the line segment (drawn with a double 
edge) between r* and c does not cross the boundary of the 
Bloch sphere. The maximal angle 9 between vectors t and r * 
is then obtained (see text) when the mentioned line segment 
is tangent to both balls. 



FIG. 8. A schematic view of the set C of vector c compatible 
with an actual state r(F) measured with maximal misalign- 
ment error e; the set C is shaded in this 2D projection. In 
our proof, we work with its enclosing sphere B(t, tX) whose 
boundary is marked with the dashed line. 



Proof. To prove C(t,s) C B(t,tX), let us take hk that 
saturate inequality (|A3I) and decompose the matrix <J? in 
Eq. flM]) as 



$ - 1 = (cose - 1) 1 3 + sine 




(A5) 



where Si = sinipi and c\ = cosifi. This allows us to 
determine the size of the enclosing ball B via the norm 
of the vector (see Figure [8} : 



d 



t = ($ - 1) t = t ($ - 1) i, t 



\t\ 



(A6) 



Since the maximal spectral radius of the matrix ^ in 
Eq. (|A5|) is \/2, the spectral radius of $ — 1 is upper- 
bounded by A [as defined in Eq. (|A4|) ]. Then the norm 
of the difference vector d is upper-bounded as follows: 



< tU> 



|*|| a <«A 



(A7) 



which shows that a ball of radius tX centered at t in- 
deed encloses all vectors c obtainable from r assuming 
bounded misalignment error e. □ 



b. Reconstructed states compatible with c£ C(t, e) 

Whenever the vector c represents a legitimate Bloch 
vector, the reconstructed state p follows immediately 
from Eq. ©: 



|c|| < 1 



P = 



(A8) 



This is true even if the state is reconstructed by the MLE 
technique. Whenever ||c|| > I, Eq. ((SJ) fails, but a physi- 
cal state can still be reconstructed using the MLE tech- 
nique by finding a quantum state p with Bloch vector r 
that maximizes the likelihood function £(r)Q where r is 
constrained by ||f|| < 1. 

By computing the Hessian of log £(r), one can check 
that C is strictly concave in r, with an (unconstrained) 
maximum at r = c. Therefore, for ||c|| > I, the solution 
r * that maximizes C must lie on the boundary on the 
Bloch sphere. In particular the line segment joining r * 
and c cannot cross the Bloch sphere; otherwise, it would 
contradict the strict concavity of C. We can thus restrict 
our attention to r that lies on the surface of a spherical 
cap delimited by tangents of the Bloch sphere passing 
through c. In Figure we plot the vector r that has the 



This can be achieved, for example, by using the iterative algo- 
rithm described in Refs. [20I Esl |33I . 
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maximal angular deviation from t while being perpendic- 
ular to the tangential plane containing c. Note, however, 
that depending on the actual functional form of C, the 
state reconstructed from MLE may have an angular de- 
viation that is less than 9. 

What is the maximal 9 allowed? Since the outcome 
vector c is contained within B(t,tX), the angle 9 is max- 
imal when the line passing through c and f is tangent 
to both the Bloch sphere and B(t,tX) (as shown in Fig- 
ure \§§. Standard trigonometry then gives: 



1-iA 



COSW = 



r-t 1- 
~ ~ t 



tx 



(A9) 



c. Worst-case fidelity 



We now recall from Ref. |34| that the fidelity function 
for single qubit states can be written as 



T(r(t),p(r)) 



1 r 
2 



1 + f • r + v/(l-t 2 )(l-r 2 )J , 

(A10) 

which is concave for any given r, and has maximal value 
1 when p = t. 

To compute the worst-case fidelity F(t,r), we need to 
consider two separate cases. For ||c|| > 1, the state p is 
pure and using Eq. (|A9j) . we get: 

t- 1 + r ■ t tX X , K . , 

2 - 2 ~ 2 v ; 

For ||c|| < 1, the reconstructed Bloch vector f = c lies 
within the Bloch sphere. To compare the minimal fidelity 
attainable in this case with Eq. (JTJ]) , we shall consider the 
intersection of B(t , Xt) and the Bloch sphere — a convex 
set which we shall denote by B'. Note that B' is still a 
superset of the set of outcome vector c compatible with 
t and bounded misalignment error e, hence, 



min J- It, c) > min J-(t, d 



(A12) 



By the concavity of £, the right-hand-side of Eq. (|A15|) is 
attained at the boundary of B' . Here, we can distinguish 
two subcases, namely, the minimizing d € B' corresponds 
to a (1) pure state or (2) mixed state, cf. Figure [Qj 

Now, we remind that the fidelity function Eq. (|A10l) 
depends only on the inner product between the vectors 
as well as their magnitude. Without loss of generality, 
we can thus write these vectors in the 2-dimensional sub- 
space spanned by t and d . For example, in the first case, 
we may write 

t = (t, 0) and d — (cos a, sin a) , a S [— a c , a c ], 

(A13) 
whilst in the second case, we may write 

f=(t,0) and d = (t- Xtcos9,Xtsm9) , 9 € [-6 C ,6 C ], 

' (A14) 



where a c = acos 



3 (±±%AV) and c = acos (£^^). 



Minimizing the fidelity for these two subcases, one finds 
that 



min F{t,d) > 1 . 



(A15) 



Likewise, in the scenario where the ball B(t, tX) in en- 
tirely contained within the Bloch sphere, i.e., when 
B(t,tX) = B' ', one can apply a parametrization similar 
to Eq. (|A14[) to show that the minimal fidelity also sat- 
isfies Eq. (|A15[) . All in all, we thus see that the minimal 
fidelity when t is a mixed state is always greater than 
or equal to worst-case fidelity for singlc-qubit pure state, 
i.e., f(e) = 1 - A/2. Thus Eq. (O is a valid bound on 
the minimal fidelity for arbitrary single-qubit state. 



2. The MLE reconstruction of a product state 
remains product 

It can be shown that the product nature of a multipar- 
tite product quantum state r = ®i=i T is preserved 
during the MLE reconstruction procedure [33[ even if 
some local errors — such as the systematic errors en- 
visaged in the main text — incurred in the description of 
the positive-operator-valued-measure (POVM) elements. 

Specifically, let us denote by r the actual state that 
undergoes the state tomography experiment, and let 
P{a\, a,2\ki, fe) be the conditional probability of obtain- 
ing measurement outcomes 01,02 for the choice of mea- 



(i) 



r(2) 



surement settings ki,k 2 , and let II„ ; k , 11^ k ^ be the ac- 
tual local POVM element used to generate these mea- 
surement statistics, i.e.: 



P(ai, <X2\k\, k%) = tr 



n 



(i) 

ai,fci 



n 



(2) 



(A16) 



We prove below the claimed proposition for the bipar- 
tite scenario. Its generalization to the n-partite scenario 
is evident. 

Proposition 2. For measurement statistics gathered by 



(i) 



r(2) 



performing local measurements Yl a fc <£>II^ k on a bipar- 
tite product quantum state t = t^ 1 ' (& r' 2 ', any al gori thm 
maximizing the likelihood function given in Ref. fSS] re- 
constructs a product multipartite state, even if the algo- 
rithm employs a different set of POVM, say, < 5^ k > ^ 

{ n ^ ,kj } f° r some 3 and k J ■ 

Proof. First, let us recall that the quantum state p max- 
imizing the likelihood function given in Ref. [331 ] satisfies 
the following equation: 



p = Rp, 



(A17) 



where the operator R encodes the maximization prob- 
lem. Our proof is valid independently of the particular 
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technique used to solve the maximum likelihood prob- 
lem encoded in Eq. (IA17J). f or example, the iterative [33[ 
and diluted algorithms [32|. The operator R is derived 
in Ref. [33l as: 



R 



E 



P(ai,a 2 \k 1 ,k2) ~(i) 



a 1 a 2 k 1 k 2 tr [/ , Pm,ti ®^a 2 ,k 



a± ,fci 



n 



(2) 

02,^2 " 



We now introduce the ansatz p 



(A18) 
r* 1 ' ® p( 2 \ and sub- 
stitute it into Eq. (|A17[) . First, the product structure 
P(ai, a 2 \k\, k 2 ) = P(ai|fci)P(a2|fc2) follows by replacing 
t = t' 1 ' (g) T^ 2 ) in Eq. (IA16|) . Then, we exhibit the prod- 
uct structure of R by introducing p into Eq. (|A18[) . giving 
R = RW (g> RW with: 



^tr[ P wn« fcj ] a ^ 

We may now rewrite Eq. (|A17[) as: 



p^®p^ 



RV®RW)(pM®pW 



(A19) 



(A20) 



and see immediately that original equation can be de- 
composed as analogous equations for the individual sub- 
systems. Since the solution to Eq. (|A18|) is unique |32J . 
we thus sec that the resulting reconstructed state must 



be 



P 



in 



p(2)_ 



/here p^> is the solution of the 



the single-qubit MLE equation p® = iJ^pO'), □ 



3. Tomography of two-qubit pure states 

a. Paramterization for low entanglement 
(concurrence sin 2a < 0.56) 

For a in this domain (0 < a < 0.29), our numerical 
results can be very well approximatcd_| by considering: 



\i> 



ia|V+>|V' 8 + >+sina|Vr>|V7>. (A21) 



in conjunction with the actual measurement bases given 
in Eq. ([3]) for both qubits, where \tjx~ ) were defined in 
Eq. @. Note that this paramctrization, in particular, 
recovers the optimal solution found for a = 0. 



the aforementioned domain (cf. Figure [4}. To this 
end, let us first rewrite 1^) as a density matrix, i.e., 
T = iV-'aX'^al an d remind that it can be decomposed in 
the basis of Pauli matrices: 



\( 



t = - l(g>l+t 1 -a<Sil + l(g)t2-v- 



-^2i T i3 a i 



(A22) 



where ti t 2 are the Bloch vectors of the reduced density 
matrices, and T is a 3 x 3 matrix that is responsible for 
the correlations between the two qubits. 

Note that for small a, t is weakly entangled and is close 
to a product state in the following sense: from Eq. (IA21[) 



and Eq. ©, if we decompose TV 



as T = t x t\ 



T, we 



see that tj = cos 2a s. Then the product term t\ t% has 
spectral norm cos 2 2a, whereas the spectral norm of T is 
||T||2 = sin 2a. Rewriting r as: 



[l + ti-a) (l + h-a) Y.i 1 Tij<Ji®<J 3 



(A23) 



it becomes clear that the contribution of t when com- 
puting the fidelity is mainly due to the first term in the 
sum. If we now approximate r by keeping only the first 
(product) term in the sum, we approximate the state 
reconstruction by solving a set of linear equations analo- 
gous to that given in Eq. ([5]) . The reconstructed state p is 
also a product, as proven in Appendix IA 21 and in partic- 
ular the reconstructed Bloch vector is (cos e + \/2 sine) tj, 
which has a norm proportional to cos 2a. It thus follows 
that the susceptibility is more pronounced for smaller a. 
We can also understand this more formally by analyz- 
ing the loss of fidelity for small e = ygj : 



L(a) = 1 - T(p, t) = 1 - tr (p r) 



(A24) 



on the optimal p and pure state r obtained from our 
numerical analysis Qj. Specifically, we write r = t s + f 
with t s = tv 2 T (g) triT, and similarly p = p s + p and 
decompose L(a) as: 

L(a) = tr[p (p - t)] = tr[p s (p s - t s )] + terms in p, f. 

(A25) 
The quantity L(a), as well as the two terms in the right- 
hand-side of Eq. (|A25[) are plotted individually as a func- 
tion of a in Figure 1101 Clearly, from the figure, we can 
see that the major (w 89%) contribution to L(a) comes 
from the marginal term tr[p s (p s — r s )]. 



b. Robustness against misalignment error for low 
entanglement (concurrence sin 2a < 0.56,) 

Here, we provide some intuition on the observation 
that pure two-qubit states become increasingly robust 
against the kind of systematic error that we consider as 
entanglement (parameterized by sin 2a) increases within 



c. Paramctrization for high entanglement 
(concurrence sin 2a > 0.87,) 



In this domain of a, the correlation term in Eq. (IA22I) 
becomes the dominating term in r and the paramctriza- 
tion given in Appendix IA3al no longer serves as a good 



This gives a 0.1% relative error in terms of the loss of fidelity. 



The loss of fidelity L(a) is proportional to the susceptibility S(a) 
defined in Eq. d 1 51) at first order for small e : S(a) sd — i2i. 
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0.1 0.2 0.3 0.4 0.5 
Concurrence sin la 

FIG. 10. Loss of fidelity in the (numerically determined) 
worst-case scenario for low-entanglement states, for small 
e = y§o ■ The normalized loss of fidelity is plotted as a solid 
line, contribution from the marginal terms tr[p s (p s — r s )] and 
the rest are plotted as dashed lines. 



approximation. Instead, a better parametrizatiorr H to 
the optimal \ip a ) and rvi that we found in our optimiza- 
tion is: 

\tp a ) = cosa \i})j) \ip£) + sina \ipg) \ipj) , (A26) 

where \ip s ) are eigenvectors of 



<?6 



sin 6 

7f 



(&x + cr z ) + COS 9 <J y 



with ±1 eigenvalues and 6 « 0.9961. The phase of \ip^) 
is such that (Oli/^) = e l ^c + and (Ol^) = c~, where c^ 
are some real numbers and <fi ~ 0.4980. 

The state \ipa) is then measured along the actual mea- 
surements axes 



,<■'> 



s e c 7 I , n 2 



C?) 



V2 






(A28) 
where s e = sine, c 7 = cos 7, s 7 = sin 7 and 7 ra 2.7946. 



4. Tomography of two-qubit pure states with 
correlated misalignment error 

The misalignment systematic errors considered in the 
main text are local in the sense that the measurement 
settings on the second party deviate from the ideal ones 

(2) (2) 

such that rh k is replaced by n k ; this deviation does 
not depend on the measurement being done on the first 
party. In contrast, let us consider now the case where 



the misalignment errors arc correlated between different 
parties. This happens, e.g. in some ion traps where mea- 
surements on one physical system also changes the state 



1.0 


OO 




Correlated * 
/ 

errors / 


1 




/ 
/ 
/ 


0.5 




Local errors 

(from main text) 



IT 



0.2 0.4 0.6 0.8 1 

Concurrence sin 2a 

FIG. 11. Comparison of the loss of fidelity, as defined in 
Eq. ()A24|) and for e — -^r, in the (numerically determined) 
worst-case scenario with correlated and local systematic er- 
rors. This graph is comparable at first order to the one in 
Figure [4] 



of a neighboring system; or more commonly when pairs 
of settings (fci, k^) are measured sequentially by realign- 
ing the measurement devices for each pair[3|. Then the 



intended fco-th measurement 



(2) 



may deviate in a dif- 



ferent way when measuring the pair (fci,fc 2 ) or (fc^,fc 2 ). 
(A27) We thus replace the ideal measurement settings rh k per- 



c k 1 ,k 2 



formed on the j'-th party by the actual settings 
which now depends on k\ and ki- 

We now compare the loss of fidelity in the local and 
correlated cases, fixing the maximal misalignment error 
at 1 degree (e = t?q rad) and numerically determine 
the worst-case fidelity in both cases. Numerically, we 
observe that both correlated and local systematic errors 
give the same fidelity drop when r is a product state: 
L (a = 0) ~ 0.025. Our results are shown in Figure [TTJ 

For small a, we have seen in Appendix IA 3 al that the 

state r has negligible correlated content T . In the 

2 
correlated scenario, the major contribution to the loss of 

fidelity turns out to come also from the marginal part t s . 

Moreover, the correlated errors seem to have a stronger 

effect on the correlated content, whose contribution to 

the loss of fidelity is negligible for small a. In this regime, 

we thus observe the same behavior for the two scenarios. 

When the correlated content of r starts to dominate 



T 



> 



1 —■ ' -*i 1 I 

El E 2 



which 



the marginal content (e.g., when 

takes place for concurrence sin 2a > 0.618), the effect 
of correlated systematic errors becomes dominant. In 
fact, the maximally entangled state (sin 2a = 1) has the 
maximal loss of fidelity in this scenario. 



This gives a 0.2% relative error in terms of the loss of fidelity. 
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